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ABSTRACT 


The two-electron hond-orhital model of tetrahedr ally -coordinated solids 
is generalized and its application extended. All intrahond matrix elements 
entering the formalism are now explicitly retained, including the direct over- 
lap S between the atiion and cation hybrid wavef unctions . Complete analytic 
results are obtained for the six two-electron eigenvalues and eigenstates of 
the anion-cation bond in terms of S, one-electron parameters V z and , and 
two-electron correlation parameters , Vj- and Vj, • Refined formulas for the 
dielectric constant and the nuclear exchange and pseudodipolar coefficients, 
as well as new expressions for the valence electron density, polarity of the 
bond and the c.ohesive energy, are then derived. A scheme for evaluating the 
basic parameters of the model is established, in which \T 2 is fit to the 
optical-absorption peak of group-IV elements in the manner of Harrison and 
Ciraci and the remaining quantities are calculated using Hartree-Fock free- 
atom wavefunctions and term values. For the twenty group-IV and III-V semi- 
conductors, we find V^- /v - Vj, '’- D but l /_ » leading to significant 

correlation effects in most properties. The theory gives a good account of 
the experimentally observed trends in all properties considered and approximate 
quantitative agreement is achieved for the pseudodipolar coefficient. Good 
agreement is also obtained for the E 2 optical-absorption peak, the die- 
lectric constant, the nuclear exchange coefficient and the cohesive energy of 
the binary compounds by scaling to experiment for the group-IV elements. Our 
calculations on the cohesive energy suggest that the intrabond overlap energy. 





discarded by Harrison and Ciraci, is an essential source of positive cohesion 
and probably rules out any major role by the interbond van der Waals inter- 
action suggested by them. The valence electron density is found to be dominated 
by the polari by and the shape of the hybrids. The preliminary indication 

is that the long-range tails of the free -a tom _S and p wavefunctions must be 
contracted to account for the observed bond density in S,' . 



I . INTRODUCTION 

The description of the valence hands of tetrahedr ally-coordinated solids 

in terms of s ( » ’ hybrid wave functions has been considered at various times by 

a number of different workers. In recent papers, Harrison'*' and Harrison and 
2 

Ciraci (hereafter referred to as HC) have presented a unified and greatly ex- 
tended version of such a theory, which they call the bond-orbital model. One 
of the major contributions of their work was the direct treatment of a wide 

range of physical properties in addition to the band structure. In this model 

3 

the band structure becomes a separate question and a large number of properties 

can be calculated with only a knowledge of the local properties of the anion- 

cation bond. Formally, this bond is equivalent to a two-electron diatomic 

1 

molecule. Harrison and HC, as well as previous solid-state workers, have 

treated this molecule in the usual one-electron, molecular-orbital approximation. 

As a step toward providing a more complete treatment of the bond, we introduced 

in the first paper ^ of this series (hereafter referred to as Paper I) a direct 

two-electron formalism. Our two-electron bond-orbital model was developed as 

5 

an extension of the method of Falicov and Harris for treating the hydrogen 
molecule. In Paper I we considered the simplest special cases of the theory 
and applied the results to the calculation of the dielectric constant and the 
nuclear exchange and pseudodipolar coefficients. In this paper we generalize 
our two-electron bond-orbital model into a full quantitative theory and extend 
its application to several additional physical properties, attempting to assess 
the importance of two-electron correlation in each case. 
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The simplicity of any bond- orbital approach rests with three approxima- 
tions. First, the appropriate anion and cation \| >' J hybrids are assumed to form 
a complete set for the description of the bond. If the j> and p wavefunctions 
making up these hybrids are atomic-like states, this set is .technically under- 
complete, although the choice of states can be optimized. In a one-electron 
description, one then has a simple two-state eigenvalue problem and the ground- 
state or bond orbital is a symmetric linear combination of two S£ 3 hybrids. In 
a two-electron description, a six-state eingenvalue problem must be solved, but 
an exact analytic solution is still possible for the two-electron bond-orbital, 
as well as all of the excited states. In either case, the bond orbitals obtained 
are then assumed to be orthogonal to one another in the solid. This is approxi- 
mately true because the four ip 1 * hybrids sharing a common atomic site are ortho- 
gonal by construction.''" Finally, all matrix elements linking the ground state 
of one bond to the excited states of its neighbors are discarded. The only 
interbond matrix elements permitted, and the ones which give rise to the band 
structure, * are those connecting neighboring bond orbitals through the Hamil- 
tonian. Then, because the valence band is full, one can make a unitary trans- 
formation from extended Block states to the localized bond-orbitals in calculating 
both the total energy and the total valence-electron density of the solid. 

Ihese latter quantities are thus given, respectively, as Just a sum of the total 
energies and a sum of the total electron densities of the individual bonds and 
are exactly independent of the remaining interbond terms. Thus physical pro- 
perties which depend only on the total energy or electron density can be cal- 
culated entirely in terms of intrabond matrix elements. 

The work begun in Paper I is extended here in several major ways. We first 
obtain a complete solution of the two-electron problem. The direct overlap 
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matrix element between anion and cation hybrids and the two-center Coulomb ex- 
change and transfer integrals, all of which were dropped in Paper I, are now 
included without approximation. In Sec. II exact analytic results for the 
singlet and triplet eigenstates and eigenvalues are obtained. We then proceed 
to develop a full formal theory of several important physical properties. In 
addition to a refined treatment of the dielectric constant and indirect nuclear 
interactions, we consider the valence electron density, the polarity of the bond 
and the cohesive energy. In Sec. Ill a procedure for quantitatively evaluating 
the basic parameters of the theory is established and application of the formal 
results of Sec. II is made to twenty tetrahedr ally -bonded solids. 



IT. FULL THEORY INCLUDING OVERLAP 


A. Eigenvalues and Eigenstates 


When the overlap between anion and cation hydrids is retained, the second- 
quantization formalism employed in Paper I loses its simplicity and elegance. 

Thus in this paper we introduce a spacial representation. As in Paper I, it 
is possible here to proceed with a basis set of the six two-electron states 
|afa > T>, |ctc4'>, |atct>, ja4 , c4->, |atc4->, and jcta^. However, both of the physics 
and the mathematics are simplified if a basis set of three singlet and three 
triplet states is chosen from the outset. Then the six-by-six Hamiltonian 
matrix block diagonalizes into two three -by-three blocks. One of these two 
blocks, that associated with the triplets, is also diagonal with all three 
diagonal elements equal. 

We begin by defining three orthonormal triplet states in the coordinate 
representation: 


!l> " Jjfp" - $ A (T, 


1 Z> = l a ' 1 ^ 4 A ( 4 > 4 )0ln) Cl U) (1 

1 >'? = let <xt>) - ri)[<Tf ; 4 <^1 


- 4 ' 
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where 5: , ( f i ■ r > is the anti -symmetric spatial functi 


on 


J-. t«;.U '> : 7-ty. [4l,in)-Ulfi)- 

n y 


( 2 ) 


In Eqs. (l) and (2) (t“) and <p c (? ) are respectively the anion and cation 

hybrid wave functions , S is the overlap integral between them: 


S =■ J 'u (r) J 


(3) 


and ( ^ and <q are the usual one-electron spin functions. We next define 
three normalized singlet states: 


'l-> 

I r / 
IO 


(ct (rh 4 ( rb 2 A • 


(U) 


( lor Ci > t \chU>) - [4(r,-> * t tih + i (Ofc'rnj 


where is the anti-symmetric spin function 


“A ■' fs 


v I trV \) ,| i ~ ) *“ 


)j (z'l 0"| v I V J 


(5) 


The singlets are automatically orthogonal to the triplets because of the ortho 
normality of q and q , but the singlets are not mutually orthogonal to one 


another: 

£ ‘ r I s > = S* 

<4 |t> - <tl(> = Sf—r)* 




^p 4r 


( 6 ) 
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The two-electron H ami ltonian operator has the form 


H (ir, , \i ^ - Up; * 1 ^ 1 } + i,; i0 


(T) 


where (~i P is the sum of a kinetic energy operator I and an external potential 

V« t 5 


Ho ( r) = T (r > + V 0 *t O . 


( 8 ) 


The quantity X zn includes bare-ion potentials from the anion and the cation 
plus the potential associated with all other bonds in the crystal. Because 
the full Hamiltonian H has no spin-operator dependence, all matrix elements 
coupling singlet and triplet states vanish, and the secular determinant has 
the simple block structure 


• !- 


0 


0 


O H I j > ~ < 0 • > £ 
i , i ■- H r, i 


(9) 


where Ej - <C I j H 1 ;) /* - | H I ? . We may write out the matrix 

elements entering Eq. (9) explicitly in terms of familiar one- and two-center 

g 

integrals. Following HC and Paper 1, we define the one-electron expectation 
values 


nc 


rl P (r) ( r ) r 


-- a , i 


(10) 
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and the (positive) transfer or hopping integral 


\\, t , ■■ - f f|> x i i ) Hr I ' '(V 1 1 " ) ^1'* 


All of the two-center Coulomb integrals have the form 


[ 11,1 l | - j 


i-i- -in , (l2) 


)i , m , i’ 1 ’ *- A , t 


The Coulomb replusive energy of two electrons on the same ion site is 


'= j ( Hi >» | H.p. y n ':= a , 


while the corresponding energy with one electron on each site is 


1 V CX ) *• ^ i S ^ 


The exchange energy is 


I (a,c l C > a ) 


and finally the transfer energy (analogous to u ) is 


H n x 1 ( n, n I a , . c ) , n i, 
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In Paper I we neglected, in addition to the overlap matrix element S', the 
quantities J and H u . It is easily seen from Eq. (12) that 3 and M a , unlike 
U„ and K, will vanish as S P , so that the theory of Paper I becomes rigorous 
in that limit. In the real semiconductors of interest, of course, the over- 
lap is large, $ ~ 0.6, and J and H p will have magnitudes comparable to U, and 
[(• Our calculations, which are discussed in Sec. Ill, show that for homopolar 
solids like the group-IV elements 


U * ?* K > H n * 'J 


with J -v -j D,. Similar trends exist for the binary compounds, where in addition 


U* > He 

because the anion hybrid is usually less spatially extended than the cation hy- 
brid. In our calculations Eq. (l8) holds for all III-V compounds except 6 V and 
BAs where the inequalities are reversed, as discussed in Sec. III. 

In order to simplify the analysis and to make contact with HC and Paper I, 
it is also useful to define the following combinations of the above quantities: 



c„ t e,„. + K « • c 

v/ r. V.. - ~z ( + H. 'I ■>' f S ( S.K + ) 

V3 0 h i £ 'CL* &:) - i ( U*.- U c )] ■ 

Si z i' [ i (u*- m ” J (19) 

V r ° ? i [ i 
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v. fa 


V - 

I 


Because of the overall charge neutrality of the solid, its energy levels will 
only depend on absolute differences of the one- and two-center Coulomb Integrals 
defined in Eqs. (10) - (lb)» 'I'hc energies defined In Eq. (19) correctly reflect 
this fact. For example, the addition of |< to £ has the effect of screening 
the bare-ion potential arising from the cation and gives t-a, a magnitude close 
to a Hartree-Fock expectation value, as shown by Eq. (28) of Paper I. Then 
( L x + t c ) is the center of gravity of the one-electron band structure and the 
term — $ ( t. a 1.6 C ) removes exactly the dependence of V.. _ on the zero of energy. 

4 . JO 

0 <• 1 1 

The energies 'J n and V 3 are conceptually similar to the one-electron HC parameters 
which we denote as \J^" and Tiie quantities y ( ^ + h* c ) and ( i..)^ - U c ) 

are appropriate average values of the electron-electron interaction, which 
physically screens V ex ^ and reduces the magnitudes of V,j and V . More specif- 
ically, V, and \L are derivable from the one -electron Hamiltonian 


H^(r)- Mr) + A 


r ( r‘ ) ta i v') J ' v i -‘ r ' v > i r 5 J 4 r 


i ( e, 


C c ) 


(20) 


—•S' ^ i 

r - r 


by the tight-binding formulas 


v/ i ' 

V° - i f«£i ~ I fa>J 


( 21 ) 
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We emphasize, however, that \/° and. y„ f enter the present theory only as defini- 
tions, All two-electron correlation corrections to the screening are included 
exactly through the remaining parameters \{ t , \'/ and \ ! {. 'the quoril.ity V /( ' Is a 
simple gen or all*/, at ion of the correlation term cK? fined in l'!q. (Pla) of Paper J, 
while \r and \\ appear as a result of the finite values of l ', , J and H a - Inter- 
estingly, even for large values of S our calculated values of and V l always 
turn out to he negligible and the dominant correlation parameter is Y /f , as we 
assumed in Paper I. Finally, we should point out the symmetry properties of 
our new quantities. The energies V4 and Vj. involve symmetric combinations 
of £„, U u » and H » as well as the symmetric quantities V 20 , S, K and J. Conse- 
quently, y ^ are covalent energies which depend primarily on bond 
length and not on the polarity of the bond. The energies V.^ and V^, on the other 
hand, are antisymmetric in c fi , U u and bf and thus vanish identically for homo- 

V r p (! 

t- and Vfe the effect 

of two-electron correlation will be to increase the covalency and decrease the 
polarity of the bond in binary compounds, because V4. > 0 . 

The eigenvalue of Eq. (9) associated with the degenerate triplet states 
can be evaluated immediately: 



c c 


k -t { 


v/ + v/; 


1 1 - 1 1 x j m 


(22) 


The eigenvalues associated with the three singlet states may be written as 


n M r - Co, + £ c - K -Y £ H ^ r| = U , T a H 

o 

where is a solution of the determinantal equation 



(23) 



11 


-v; -v l in 
-v; <i;; 

v/ -*v; v-ii; 


;isv‘ VT '*E M -V/-‘.VV vv 'Vh\ 

; v;”- I:m “V; i.svV V'" ^.h 


(24) 


Equation (24) is, of course, just a cubic equation in and reduces to Eq. (20) 
of Paper I in the limit S = - V^° ~ 0. As in Paper I, the physical content 

of the solutions of Eq. (24) is most transparent in the two special cases where 
this equation factorizes. These cases are: (l) \T° ~~ \ r °-~ 0, which is the 

appropriate solution for the group-IV elements, and (2) < V ? “ vl - o, which 
is the limit of no two-electron correlation. Following Paper I, we denote the 
former case as the Falicov-Harris limit (hereafter referred to as the FH limit) 
and the latter as the Harrison limit. 

In the FH limit the energy eigenvalues and the splittings between them 
may be expressed as 


i= * * 2 a -k + >sv a + v; - 2s\ vi>-Vt)/o~/) 

E iv. - -t; -K + z s Vk t2V,. *-~vl - 2s( v„-y, 

Ex, a - 2 £ - k + 2 s\-j + \ ± (4V/+ v; s )t. : 

and 


E K - E„ zt - 2 4 V,(i-3sVn- ; i- ('l V/ I VV' ) 

" E £ - - Vi, ( 26 ) 

Cj - e H - 2 . V1 / U - S . 1 ) - ^(i+s ! )/ 2 -:;')t( 4 W< V) 



where we have dropped the subscripts a and c on £ and have introduced the re- 
normalized variables 

V, | V. f *> ( v/) ~ V- ) j 7 ( i -.s'" ■) 

V; v V/i-i')*- 

\ - %'/ ( J (27) 

Vy = v c "/ii-y 

v,. . Vt'/ f.-V) 

which facilitate comparisons with HG and Paper I and which will be used in the 
remainder of this paper. The renormalization factors of ( I — $ 2 ) ^ and (l — $~) 
are precisely the same as those used by HC in defining y^and respectively. 

The term added to \[^ may be viewed as the subtraction of the constant V 4 - 
from H 0 in Eq. (20). The qualitative ordering of the energy levels in Eq. (25) 
is the same as given in Fig. 1 of Paper I, with E the ground state energy and 
the triplet levels and the first excited singlet ,£^ separated by Note that 

overlap has shifted all six levels upward in energy by an amount 2 S\f x This term, 
however, will contribute only to properties which depend on the absolute position 
of the levels, such as the cohesive energy. For V(, ” 0 the triplet energies and 
are also shifted downward from the center of gravity by 2 S“ Y|_ /( l~ S 2 ) • 

This means that if one fits Ejg - Eyj_ to the principal optical absorption gap 
for a fixed V^, as we did in Paper I and as we do in Sec. Ill below, one requires 
a larger value of Y ? than with S ^ 0. . p 

In the Harrison limit the energy levels and splittings are found to be 


~ I. r, 1IL ™ c Ji “ £*\ ’ h K + 2.1 \/. t 


( 28 ) 


h 


T, tzc /“ £<\. + £ c ~ K + 2 S V z i 2 ( \ 2 \Y 


7 •’ \ 'M 
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Again all levels are shifted by the amount ZS\^ which is consistent with the 
shift of S \' in the one-electron energy levels of HC. Otherwise, overlap has 
not altered the formal structure of the eigenvalues and Eqs. (28) and (29) are 
completely equivalent to the HC results, as discussed in Paper I. 

We next consider the calculation of the eigenfunction of the Hamiltonian 
H. It is evident from Eq. (9) that the three triplet basis states are already 
eigenstates of H: 


IX? = I ' > 

(30) 

| If/ \i> 

I III > “ I S > 

The singlet eigenstates can be written as a linear combination of the three 
singlet basis states: 

I m > - £- a lw H > , r m-- 1 V (31) 

The ground state jH> is the two-electron bond orbital. Setting 


= A iH /p ; i 


( 32 ) 



l4 


where [) is the normalization constant 


L; 5 '\v/m A r M + A^I + - 2 - fj A'i^AcM ^ A n.M *' ) A<M j ^ 


(33) 


the coefficients Ay^are obtained in the usual way and may be expressed as 

; w '.= pi f n-2SV--ViM<'^) •+ sshv.-H)- s‘£m j£-V,0-i‘> 

A 1 « . i H s'E 1 - 

+ Stv-, -VO 1 ; S'/,* ■Vr'lC-l; P .* - Sflp-£2V + 0-sPi2V 3 (i-i ! P 

- E,” j[-V 2 ii-s* V + StWVt'i f (w,+ v J ;Vi-s z P -SEm ] j (34) 


A/,, 


6 m 


-• v, 1 1 ;P- C< -sv, -ViYd-s") 


+ 2$'-( VS,-Vii- & Em J 


In the general case, once the roots of the cubic equation (9) are found, they 
can be inserted into Eq. (3*0 to find the 

In the EH limit, - D — 0 for M - TT and the general solution Eq. (3*+) 
cannot be used. In this case, the proper coefficients can be shown to be 


f a -l - : I / 

^ ft f IT 


(35) 




For M, -V , YTj on "the other hand, Eq. (25) may be used directly in Eq. (>34) to 
obtain 


G Els 
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& M .KI ' 


J 


6 m 


w,, -- 4 [> ('•+ M ) t SGJ/P, 
'MM-- ' 


r x o-^" ) i- Q (i+ ) " y- ] / 


(36) 


with 


D, 


\ ■ X 


■*. ! J J 


r 4 - 


■A 


f X 1 j (i-S'M 


(3T) 


where we have defined 


x - - v-. ■ v* 


(38) 


j. 


(39) 


and 

Q v ( X ? '> 

Note that all explicit dependence on has cancelled out in this limit. 

In the Harrison limit the coefficients are most conveniently expressed in 

■^2 

terms of Harrison’s polarity and covalency parameters ’ oL and ok.., respectively 


rx 


and 


CA, 


Y.> / ( \ X 1 ;Y ) 


(Ho) 


Wt v; 

( )- y i; 


(Hi) 


One finds 


f ^ 4 IV' / 


l 


a. 


-+ j/O-s" ) 


(42) 


‘ .it J 


a 


1 1 -- 


/0~Y) 

^ 1.3 1 




P00R 
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and 

- ..u 
I a i 

^ (- M 


X :i \» 
M 


•«, | i d (• t s ) • ( £ v 1 ;; > j / \ I - r ! ) 


(U3) 


si 


. |! s; -f s ~»'J /tbs?) 

Z - ' 


■where Q- - +J for M - V and 9 ■ ~l for | v j - .11 . It is strightforward to show that 
all of the above results for the FH and Harrison limits are equivalent to those 
of Paper I in the limit 5 ’ \ r {, - 0. 

B. Physical Properties 

The full effects of two-electron correlation and overlap on a wide range 
of physical properties can be assesed 'using the eigenvalues and eigenstates 
obtained above. We shall concentrate here, however, on only a small select 
number of such properties. For the purposes of comparison with HG and Paper I, 
we derive generalized formulas for the dielectric constant e and the nuclear 
exchange and pseudodipolar coefficients | ^ and f' , respectively. As an 
extension of the method, we also consider the valence-electron density, polarity 
of the bond and the cohesive energy. 

1. Dielectric Constant 

As in Paper I, the quantities f , and I ' ^ are most conveniently calculated 
in perturbation theory. Following HC and Paper I, we apply an electric ileld 
($ in the +x direction and consider the energy shift A E induced in an isolated 
bond lying in the [ill] direction. The origin of coordinates is chosen at the 
center of the bond and, in analogy with Eqs . (3k) - (J+l) of Paper I, the shift 
in the total energy (per unit volume) of the crystal to second-order in E is just 



IT 


y 

k-r. Js 


A j 


V 


(I!-- ) 


t i 




(44) 


where the polarization is 


P, 


( ) n a { ‘'f i '•! " c 


(45) 


and. the electric susceptibility is 


% 


n y 


.1 N 


i' l IL’ , r 


l i M A H ; 

‘ v i i H 


(46) 


Here j \ *i i ^ Im /^ { 11/ is the ground state, M is the average 

valence-electron density, and (J is the vector distance from the anion to the 
cation. Note that the sum in Eq. (46) runs only over the two excited singlet 

O 

states, because U c cannot couple [ (q / to the triplet states. Using Eqs . (U) 
and (3l)» one can derive a general formula for j l(^ |(q / in terms of the 
expansion coefficients 

- ( d / {*> ) _ l i ‘ ~ ^ 4^ 1]-|) 

~' r 1 y j jVA \ A, ' ^4vt Vl'M + CQsu- ditv ] j 


-( 4.. - [\- v ) 

1( a ri> 

c - 1 

^ i- H 

■+ vl/; tv| ) i' 

y r j A y m 

/h ; ) 

L -J- >. :. ;; L 

2 2 - J V lts : ' 

[ 

+ 4,; 

f i,i T ^ ' A ' i ! 

1 M -tf r-i ) 

* ] 

_ r* r , 

v 2 ) ^ 

J C(K 

j r ( 


! ~i y ; i , L k ,\ 1 ) 7 

fc, 41 k i - 

/(H C ) j 


^ 2. S j-",, [_ ( - 1 ' 2, U ' ( ) ilj K[ -t ( ! '> '• i ■' 4 -/ r.-i ) r j j 


( 47 ) 





where we have defined 





I I X M '| ’ n i 



) 




m 


and 



In Eq. (hS) X'n is the X coordinate measured from the center n. The quantity 
[ I - (8.\ f Sv." 5 j is the f 1 defined in Eq. (l 6) of HC and Eq. (36) of Paper I. 

The renormalization factor of ( I ~ ) 2 is again such that our f is precisely 

the same as that defined hy 1IC. Note that b n cj is just the center of gravity 
of the hybrid electron density A Anew overlap term $ c f not consid- 

ered in either HC or Paper I, has also appeared. For homopolar solids <ica- P 
and S c - £ "by symmetry; for binary compounds one usually has o c > and dV x > 0 

because of the greater spatial extent of the cation hybrid. 

In the FH and Harrison limits, f* 0 . and OC are given by particularly simple 
expressions. In the former case $ ~ S - <L = 0 811(1 one finds, as in Paper 

I, that only the matrix element j T)[ V is non-zero. Thus p 0 ” 0 , as expected 

by symmetry. Using Eqs . (35) and (36), it is also straight-forward to show that 
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T ’ Vs ~ 


•1 3 V N ( d r f \( v. ) 


( 51 ) 


with 




Q (f 0 X ') V *+' y (i- 


) -f iv ( f 


/V' ) 


(52) 


where we have defined 


V,./ v. 


(53) 


Equation (5l) has the same analytic structure as obtained in HC and Paper I, and 
the factor & reduces to Eq. (46) of Paper I in the limit 5 - A ~ 0. In addition, 
fi " \ for V - 2 0 and j> > | for 2 > o and \ 2 C. A further simplification is 
also possible. In HC the quantity V^ 1 is determined by fitting £ ^ - jZ ^ 

principal optical-absorption peak. If one uses the same standard to fix V^, then 
one has, comparing Eq. (26) and Eq. (29) with V-, - 0. 


:s. Vx l ~ ~ I 


XO-Yi (i-s 5 Yj /i'i-sv) , 


(5 1 *) 


and the effect of the two-electron correlation on the HC result can be expressed 
in terms of the simpler factor 


?' t fv v; : -; - ft x y, (55) 

which depends only on V and not on either S or 2 explicitly. Our calculations 
suggest X~ I f or the group-IV elements, in which case ^ ^ !. >. In addition, 

if the quantity is determined by fitting € to experiment, then the fitted y 

■■ i;.- -■ 

will be larger than the HC value by a factor of ( f 

In the Harrison limit IF ? and <(<qj U^l {*•/ are finite and j TJ> 

vanishes. Using Eqs . (29) s (42), and (43) one finds in this case 
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V 
1 f 


v v-n 


) ( < • \ ) 4 / 


A ■ 


sU, J/6W >7 7 


(56) 


and 


c. 


4 - *!( Ti fj v 'v 


* / 

Vi. 


•I r \V7.C v, • w / 2 


(57) 


(58) 


where we have defined 

f ^ + (‘•M./fxc } _ ;.v V) “ •" )- ! _ 

Equations ( 56 )— ( 58 ) reduce to the HC results in the limit L - and a ' 0 , 
as was implicitly assumed in their work. Setting $ = & CA ~ °> one recovers Eqs . 
(48) and (51) of Paper I for and £ , respectively. 

We should mention at this point that the vanishing of the matrix element 
<\4 I V only appears to be an exact result in the FH and Harrison limits. 


However, our calculations suggest that the additive term |<6| | i W | /f£4 - ) in 

Eq. (46) makes a negligibly small contribution to X. in real materials. Thus 

, . He 

even in the general case, the correlation correction to 'X is essentially a 
multiplicative factor as in Eq. (55). 

2. Electron Density and Polarity 

The total valence-electron density in the solid is just a sum of all the 
individual bond densities. The bond electron density may be obtained from the 
two-electron ground state f 6j )> in the following way. If the spatial part of the 
singlet basis state |t> is written as /l/o(r, f 1 ) , then the electron density asso- 
ciated with the ground state is from Eq. (3l) just 


, Cl ^T a j-i . 


r> r') if;. [ r, r'V f ' 
T v j 


p ( r) -2 2 . 

’6 ) - h, 1 , h 

The factor of 2 arises because there are two electrons per bond, i.e.. 


( 59 ) 


j f (r)c(r = 1 . Also note that ^ > r ' ) 


OkIGMAE PAGE 
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is symmetric in f and. p , so that it does not matter over which, coordinate 
one integrates in Eq. (59)- Using Eq. (4), it is simple to show that 


f (r) - 

l -jV‘ 

( r \ + 7(c Unc 4ii<-r\ ' : K ch 

(60) 

with 





- a 4d 

f «4' / < / + $“ ) 2 'Zyh- a ^’ 

(61) 

■«r ■=• 


+■ ^k>C\ /( ! '4- ) + -- • |,'ri v UtT 

( 62 ) 

and 





Hp - n 4u Ju, <• 


, I- S«46 /(» ■) 1 vlhu 1 Mu, 


(63) 


where we have noted that { i '■ and r) are real. In the FH limit one 

finds using Eq. (36) 


tr : L“ : 


1(i = 4 1 


i s ; ^ t x V x ); ( i - )/ 1 u _ 

(5 (t-$$ )i)+ sr'>.--5Z & - C 'HA - * ) + $>XQ 3 Xp 1 ' 4 Xp . 


In the Harrison limit, on the other hand, one is led from Eq. (43) to the results 


■lU I 


«.<r 


2 »■ 


= [(i - i x i< ~ J ' . , 


i-o ) 


(65) 




and 


p OOp P 4 q* 


uz * X\ iu ~ 


( 66 ) 
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For X ' 0 in Eq. (64) or p ~ and d c ” f in Eqs. ( 65 ) and ( 66 ), 
llj l(/ ~ Va\ • (H S )' ' > ns one might expect. 


We tenutively define a polarity associated with the electron density 


as 




<4 

r 


s i ( Hi - «/ ) ( >-s' . 


( 6 ?) 


From Eq. ( 65 ) one clearly has d p ~ p in the Harrison limit, so that Eq. ( 67 ) 

is consistent with the HC concept of polarity. In the limit of no overlap, 

both f; = 0 and ^(r'j^^^Oand this definition is relatively unambiguous. 

In the case of large overlap, on the other hand, the definition of o(p seems 

less compelling. Ultimately, an arbitrariness arises from the innumerable 

ways one can divide up the electron density in a periodic solid. Nevertheless, 

4 

there need be no additional uncertainty in any calculation, so long as o <(p 
enters solely as a definition. 

The concept of polarity is also intimately connected with the polarization 
of the bond P 0 . Using Eq. (47) withM = <q- in Eq. (45), one has in the general 

C 

case 

P p "" <t > (68) 

where 

(69) 


"Z p ~ Y*dj? f ~i' C ld<\ t Kc ) •- 'll 


2 . 

a c 


The quantity is the effective charge associated with the dipole P 0 . It 
would follow, therefore, that Zp is also a reasonable measure of polarity. 
However, this definition has the disadvantage of depending on the additional 
parameters Y ^ } and In HC the quantities £ c - a^and S^are set to 
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zero and ’j* is fitted to experiment through the dielectric constant Eq. (57)* 

The fitted \ thus adsorbs local-electric -fie Id corrections to £ as veil as 
corrections to the bond-orbital model itself. Consequently, it is not appro- 
priate to use their in any definition of polarity. On the other hand, one 

7 y \ 

may calculate 1 y directly and compare it to J. y ? and we do this in Sec. III. 

We finally point out that one may, as we did in Paper I, define a polarity 
for each of the states / NJ > by an obvious generalization of Eqs . (59)-(67)* It 
readily follows from Eqs . (35) and (42) thatot-p = 0 for all three singlet states 

, IX i T , (( 

in the EH limit and that oCp 0 and y p •" - <Xp in the Harrison limit. These 
results are in agreement with those obtained in Paper I. 

3. Indirect Nuclear Interactions 

To obtain generalized expressions for r e and , we again consider the 

. — ^ 

magnetic interaction between an electron with spin \ (T and a nucleus with 
spin ) . In analogy with Paper I, this interaction couples only singlet and 

triplet states and to second-order gives rise to an energy shift out of the ground 
state of 


4E = - 1 - > ' <^|1<a Ih XmI Ua / d > , 

^ *- t M " j - ;I , Id 

where 

W A ( i ,, t. )■ - Z. JL Viirii-d n) 

c-i,a. n-a'jc 



(70) 


(71) 


( 72 ) 
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In Eq. (72) J ( p is the Bohr magneton, Y t1 is the gyromagnetic ratio and P is 


the electron coordinate measured relative to the nucleus u . The triplet states. 


as given by Eq. (l), have the form and the sum over y\ in Eq. (70) 


-A ’ ’ ' ' " M 

can he immediately accomplished by using the sum rule 

1 < Z A | <r ( i 1 1 h. M >< 2 m I r (0 1 Z A > * 


it i ; p* 


(73) 


M-- I . U> !£ 




(-0 I 

Then using Eqs . (2) and (31) and noting that ^ (P 1 P) is antisymmetric in p, 
and ip , one may express ^ £ in terms of the expansion coefficients 0. > ( . ■ : 


r - 1 V 

A r ^ ~ 7- > 

fc ° n ft, C. 




■n L\ >\ \ ) 




(7 1 *) 


where •' 1 and 0 - - / and we have set 


r ,, « f\. J < 1 1 A.. l«h A > j L i a > d 1 " ) / :: ( ,/ * , / A„ 1 -p , 


t!A„ip'-) ' i-t I A„i'i> c >] ;/(/ 5-’ ) 


(75) 


b, " V ( ~ v,> i 'I ■) 


(76) 


b, l •'> ( “p t . * "b "bti ) 


and 


■e: 


D > 


(7T) 


p r*~ » 

b n - /•£, 


Equation ( 7 U) is a generalisation of Eq. (62) of Paper I. In addition to S , 




four new overlap matrix elements have appeared in A £ , namely, ^ ] /\^ j -,p > 
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<kj Ac { t '> » \<t ; he \ > , and <<P C I Aai ( l c > • It is not hard to show that all 


matrix elements of A, u have the form 


<1,1 Ml > 


/“p 1 *VL i 


II 

1 1 n" 


1 * 


r, it 


vi, r> , >1 


(78) 


with 


<£-* 

IT 


A A 

3 fac ) fl c 


- 1 


(79) 


and where we have defined 


6 TT < 4v l ^ < P) I > 


( 80 ) 


and 




(81) 


In Eq. (79) r A c i s a unit vector directed from the anion to the cation and in 
Eqs. (80) and (8l) the electron coordinates p >v / are to he measured from the 
center V . Also, the right-hand side of Eq. (8l) is written for the Z axis in 

. / t 

^ A V v^> 11 

the direction of . More explicit expressions for and ,, in terms of 

the S and J components of the hybrids are given in Sec. Ill and the 


Appendix. As expected, the non-overlap terms ^ and ^ . immediately reduce 

6 * * 

to Eqs. (64) and (65) of Paper I. Furthermore, the tensor Q n clearly has the 
same form as Eq. (78) with ^ " ^replaced by J ^ 




" ] 
J flC i 


1- r ’ ) Yi 


(82) 
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and $ , replaced by 

in - >>»' l ( I 1 '• j,' ) X, -J ^ ~ 


( 83 ) 


* > * > 


** V 


Writing out in terms of and )P M and noting that ir * TT = £ + "TT , 

we arrive at the final results 


i,P 


where 




, / 


* . T 9,?,: x„-( p e - . * + ^ j-i„, ; 

Hj ))*-'- 4|<1 


(8k) 


II 


• ~ / { r, it )’► rv L tV i»i»' + * • ■' / 


( 85 ) 


and 


r, 


)i •' 

M 




; V t. vv [ iu C ,- 4 ’ 4 & y .--*A' 


( 86 ) 


Equations ( 84 )~( 86 ) have the same analytic structure as the corresponding results 
obtained in Paper I, with all effects of two-electron correlation and overlap 
being absorbed into the quantities ^ and E 0 • Again we identify Q, = Q. 

and -• jlj as the nuclear exchange and pseudodipolar coefficients, respectively. 
We normally expect and )p n and consequently and to be positive 

quantities, as is observed experimentally, provided 


0 £ h >/b 2 ± $ 


( 87 ) 
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and 

c - ^ ' A-‘- 

c > ■' 

for »'/ n . The lower limit in Eq. ( 87 ) is exact for homopolar solids, 
where - We find Eqs. ( 87 ) and (88) obeyed for all group-IV and III-V 

binary compounds. 

In the FH limit, Eqs. ( 82 ) and (83) simplify to 
The coupling constant £ * can be shown to be in this case 


( 88 ) 


(89) 


E. 


4-V, 


< ) 


(90) 


where 


f , ;4('Q+ X-+s) 2 l | -S^/[0(R+)Ot4-S(a-y)rt^- XOtS 2 ) (91) 

which reduces to Eq. (72) of Paper I in the limit Also, J ? I for £ >p 

and/or X > 0 , and ^ =? | for )( = 2-0 . If V z and 1/^ are related by 

Eq. (5*0, then the HC value of £ g 1 is enhanced by the factor 


7 '= v,"' 7 / W - il-;z + X(i- 3 s 2 )+Q 0 -^)]j / fs'), 


(92) 


which unlike b' is a function of £ and i as well as Y . In Fig. 1 we have 
plotted Z 1 VS X for 2 - P and several values of ^ . Clearly, 7 > is a sensitive 
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function in the region of physical interest: 5 /v - 0. k , X •' v - /, 0 

sensitivity results primarily from the energy denominator in Eq. (77)» 
which can vanish for sufficiently large X . 

In tho Harrison limit 

v - a„ i t * s (otf./o u )j i;; - <V*v>)2«: 


This 


2- ( j aUc ^ /(' S ) > 


(93) 


with a similar expression for , and 


- 1 


_L 

A- 


Vi/c v; • vy ) % . 


(94) 


E. 

Equation (94) agrees with Eq. (73) of Paper I, "but note that for J[ u ‘t . J| ^ --J! c - o 

and. ^ AO . We do not recover our previous formulas for and p ^ h " * 

Contrary to our original expectations, all terms involving ^ cannot simply he 
absorbed into Va. and . Both ft 1 '” and P " are multiplied instead 

by a factor of [ |~S ) in this limit. 

4. Cohesive Energy 

To obtain the cohesive energy, one needs both the binding energy of the 
(two) bond electrons in an anion-cation pair and the binding energy of the (eight) 
valence electrons in separated atoms. The latter can be written as 


z(£^ ■+ ) f (i+*?) ffP t (2 •''2)£p7 A?=c,.,u 


G *P 


E 


. f 


(95) 


at, 




2 c 1,1 £ . ,e t 5f (f 

^ C scv c £<- • 5 O pa 


At 


where A 


s *7 ^ 


for group IV, - | for III-V, - 2 for II-VI and — 5 for 
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I-VTT elements. Also -In Eq. (95) 4. t , and are defined an the binding energies 

of the v ninl )' e lor Irons of atom H . We use the supersirr i pt, ill' 1 here bemuse 

H> m 

in practice I t L :v convenient, to take t,. u and to be liar tree- foe k. 1'ree- 

atom eigenvalues, although in principle this need not he assumed. The corres- 
ponding binding energy of two bond electrons is just our ground-state energy !-~.~ 
plus the electrostatic energy of the compensating nuclei: 


r 


fcw'ft 


= tr, •+ 


e '7d 


HP, HP. , r 1 

- r e, + ^ 


c 

4 


(96) 


where we have defined the hybrid energies 
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;u 


, , Hl ^ 

-* c. pt' ’ 
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(97) 


(98) 


(99) 


and the electrostatic overlap energy 
with 

r iff 

* t y) ■-*; - e H > * a-, c . 

1 r- 

Both Harrison and HC discarded in their treatments of the cohesive energy. 

However, for the case of large overlap this is not justified, since we expect a 
to be small and negative and K to be significantly larger than e /j( • 

Using Eqs. (95) and (96), the cohesive energy (per atom pair) is just 


E 


Cxr t 


' Z rt.X'r»H.u 

— - 4- E 4 
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f - E ~ t i uin _ 


( 100 ) 
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where £; , ,■> iu the promotion energy defined by Harrison^ 

I 

c k i » i v" -^) vr v 
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and f i 10 ., a is the transfer energy 
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( 101) 


(102) 


(103) 


Both Harrison and HC took the transfer energy to be 2 a / V_, . Prom Eqs. (l9), 
(27) and (99 )» however, one can see that this is only formally valid here if 
1'(,v Ut * A ~ J ‘ s Sc and S ^ C • In the FH limit with ^ 0 and 

E'p u = B V, " » one has 

r a I ( A V " + y ^ j " - ~ ■" V, - V4 "F V(, j - f El * - K V, 11 

" 4 L(4V ' ^ ‘ J 1 • (io4) 


In the Harrison limit, on the other hand, one obtains 


£' j ' ~ 4- [7 ( V7' t \V Ya ~ 2 ^ J ' ^ ^ ** £^ A1 . 


(105) 


For E (, ( •' r and £ 4 - 2 A z? , this result agrees with Eq. (34) of HC. 
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Physic alLy, we may interpret the various contributions to E ^ as 

follows. The promo bi on energy is the energy required to promote eac.lt valence 

» 

electron form i t.s nUvuilc .!» or l’ orbital to an hybrid orbital. The 

transfer energy Is tile energy required to transfer & i. electrons from the 
anion to the cation at infinite separation, so that the four hybrid orbitals 
on each atom contain one electron. The atoms are then brought together and 
bond orbitals are formed from the overlapping hybrids. Occupation of the bond 
orbitals returns an energy - + E 0 f) P er l° n d - Equation (100) thus correctly 

includes all intrab ond contributions to the cohesive energy. To be sure, there 
are neglected interb ond contributions, as there are to all properties one cal- 
culates with the bond-orbital model, but we certainly expect these to be of 
lesser importance . In HO the possible importance of a van der Waals inter- 
action between the bonds was argued, but this was introduced in an ad hoc 
fashion to explain, their calculated negative values of Berk . In contrast, 
our calculations suggest that there is no fundamental difficulty in under- 
standing the cohesive energy of tetrahedr ally -bonded solids in terms of Eq. 
(100) alone, although the precise results one obtains are somewhat sensitive 
to the details of the calculation. 


Q r 



III. QUANTITATIVE APPLICATIONS 
A. Parameters of the Two-electron Model 

We now consider a detailed application of the formalism described above 
to real semiconductors. Our analysis will be considerably more complete than 
that given in Paper I, as we seek to understand quantitatively both the relative 
importance of the various parameters which enter the two-electron theory and 
the effect of overlap and correlation on the predictions of the one-electron 
theory. We must first establish a procedure for evaluating the basic para- 
meters of the two-electron model. We are guided here both by the experience 
of HC and the fundamental limitations of a bond-orbital approach. One must re- 
cognize the approximate nature of any bond-orbital theory, so that the difficult 
task of a complete first-principles analysis is not really warranted. On the 
other hand, one can reasonably expect the theory to produce correct orders of 
magnitude and significant trends for a wide range of physical properties. Thus 
it seems desirable to fit only enough parameters to provide a proper scaling 
to experiment. 

Of the various parameters which enter our two-electron model, one may 
separate out those which depend self-consistently on the full Hamiltonian of 
the crystal, and consequently are difficult to calculate accurately, and those 
which depend only on the specification of the hybrid wave functions . In the 
former category are the quantities V, g , a £ a, and A £ c . Both and 
enter the covalent energy ]/~ z , while the polar energy depends upon the 


- 32 - 
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difference a£ c - & £ a . Following HC, we fit V 2 to the principal optical- 
absorption peak E, for each of the group IV elements. From Eq. (54) with \V- r 




0 



l i ’ - V 


.1 /2 

( 

J s 


(106) 


which is easily evaluated once £ andV^ are specified. It is then assumed that 
\ r z is constant for an isoelectronic sequence (e.g., Ge, GaAs, etc.) and that 
values for skew compounds are given by the appropriate geometric means of the 
group-IV values. The quantity was fitted by HC to the experimental di- 
electric constant of each binary compound through Eq. (57)* This is less 
convenient in the present case because V3 enters our theory through the rather 
complicated cubic equation (24). We have chosen instead to make the assumption 

,4 £ c " - 0 (107) 

and then to calculate as 

V, • i [ ( C - 'i (M»- UiV] / (108) 

Since i A t ,-, / £ „ ( is small, one may reasonable expect Eq. (108) to be a good 

approximation. From Eq. (l8) it is clear that the term Uc) will reduce 

the magnitude of V3 , so that normally V3 < 2 ( f/* - ^ ) > as was found 

empirically by HC. Interestingly, the two exceptions to Eq. (l8), namely Bp 

and BAs , are the compounds for which HC inferred negative values of V^from 

■ , ><> nr 

their fitting procedure. From Table I it can be seen that the quantity -j; ( - £ a ) 

does become much smaller for BP and BAs than for the rest of the III-V compounds. 



However, this effect is offset "by the change in sign of the term -j ( l|„ M , ) , 
no that norma 1 punitive values of \/ s result from Eq . ( 1 OB ) . 

The remaining pm* functors V/, , V, ul \J h ure functions of the one- 
center integrals 1 ^| IV and the two-center integrals $ , K , Hn. and J 
All of these integrals depend only on the hybrid wavefunctions ^ ( r ) , so 

that we also choose to calculate these quantities directly. In the spirit of 
the bond-orbital model, we construct our hybrids from the appropriate Hartree- 
Fock free-atom J. and _jt wavefunctions. We shall not entertain the very diffi- 
cult question of exactly which atomic orbitals constitute an optimum basis set, 
but rather we treat our choice as an additional assumption to be tested. Choosing 

A A 

r « 0 £ » one con write 

tir\) - - [ R lrt (nt + 3 Gn. R f .(r,J ~ ] , (no) 

— ♦ 

where F a is the electron coordinate measured relative to the center to and 

K Sa ( r n ) and Kp^vCrn) are assumed to be positive as f a 00 . Free-atom 

Hartree-Fock S and £ radial wavefunctions, r a R i>v ( r a ) and r K Q ) » 

as well as the corresponding term values £g F and , have been calculated 

7 

and tabulated by Mann for the entire Periodic Table. The use of such atomic 
tables is convenient, but it does restrict one to consideration of compounds 
formed out of group III, IV and V elements, since the states of elements in 
groups I and II are unoccupied and are not usually calculated. The details of 
evaluating "U a » 5 » K » and J in terms of the hybrids defined by 

Eq. (110.) are discussed in the Appendix. Briefly, U -n- can be written, as a sum 
of the f and (7 integrals defined and calculated by Mann , as shown in Eq. ( 77 ) 
of Paper I. As a test of our numerical procedures, we have also evaluated these 
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p and G integrals directly from Mann's wave functions . Our calculated values 

of U H are listed in Table I and are to be compared with those given in Table 
I of Paper I, which were inferred from Mann's tables. The agreement is better 
than 0.1%. The two-center integrals S and H )x , , on the other hand, are most 
easily calculated by expanding the S£* hybrid from one site in terms of 
spherical harmonics centered on the second site. Similarly K can be evaluated 
by expanding the Coulomb potential arising from one hybrid electron density 
about the second site. In both cases, this procedure leads to a finite series 
of one -dimensional integrals . This same type of expansion method applied to the 
exchange integral J , however, leads to an infinite series of terms, which must 
be truncated. Tests of our procedures in the case of hydrogen IS orbitals, where 

g 

exact results are known , suggest that we calculate 6 , and K to an 

accuracy of better than 1%, but that we may underestimate J~ by a few percent. 

Values of 3 > K and Vy that we have computed for the twenty group-IV 
elements and group-IV and III-V binary compounds are listed in Table II, to- 
gether with values of ]/ 2 obtained from Eq. (106) and Vj obtained from Eq. (108) 
The two-center Coulomb integral K decreases with increasing bond length, as 
expected. The overlap integral 5 , on the other hand, tends to be constant, 
although AIN, GaN and InN have somewhat smaller values than the rest. For the 
solids listed in Table II, one has approximately K ~ ()+$> fi/d ~ /. 4 £ /ct 
The constancy of S and the simple dependence of K on bond length give V^. its 
expected covalent behavior. The large magnitude of S » however, alters the 
relationship of AT Z and ~V Z C anticipated in Paper I. For $ = o Eq. (106) 
demands that -c V 2 Wc \ but in our case g 2 > and we calculate 

Vz > Ya . For the group-IV elements X ~ 2. VqJ Vz - 7 7 for C 

= /. 0# for Si, = MO for Ge, and = /.IS for Sn. ■ 
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As mentioned in Sec. II, our calculated values of and V& always turn 
out to be small in comparison with V,-, and V/> » respectively. Specifically, 

we f i tul 

0,11 > £ - 0, 0 2 C V (111) 

I Vjr/Vjl < 0 . 01 ^ 

and 

o. 27 > Vfc > - D.02 £ V (112) 

I V t / V, I ^ o- > 1 

for all twenty solids considered here. In view of the large overlap of the £ 
hybrids, the uniform smallness of Vf and V4 is indeed remarkable, and we have 
found no simple explanation for these results. The magnitudes of Vr and 
are actually comparable to the numerical uncertainty in these quantities. Far 
the sake of conceptual simplicity we shall set = l/t ~ 0 in our subsequent 

analysis . Note that with fit to experiment and V? and set to zero, 

values for the integrals j-Jn, and T are no longer required. 

Having established values for the basic parameters of our two-electron model, 
we may quantitatively solve the cubic equation (2k) for the singlet-state eigen- 
values. We have done this and the values of , defined by Eq. (23), are given 

in Table III, together with the triplet-state eigenvalues 

Ey = i sV, - ■is i y+7.(c : s i > ' 113 < 

for the twenty solids under consideration. The corresponding eigenstate expan- 
sion coefficients Cl 4M , CLm and C( m , as calculated from Eqs . (32)-(3^) > are 
listed in Table IV. Also, in Table III we compare our theoretical values for 
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the principal optical -ah sorption peak E 2 = Ejj. - with loth the HC pre- 
dictions, ~ 1 [( V/° ) + ( ) J , and the experimental E^^and 

9 r # 

suh-peaks as given hy Phillips. Because for the group-IV elements rep- 
resents a fitting parameter, only the results for the binary compounds offer 

f- He 

a test of the theory. Both our b z and for these latter materials agree 

with the available experimental data to within 12$. Our &2 tend to agree 

r r 

best with fc: z g and fc 2 with , but the differences do not appear to be 

significant. It is important to stress, however, that our , in contrast 

_ tie 

to , include no direct experimental data for the binary compounds. 

B. Computed Properties 

With the information given in Tables I-IV, we may systematically study the 
physical properties discussed in Sec. IX. We begin with the dielectric constant, 
for which we need the additional quantities "T, )T C » and 5^ . All of these 
parameters are readily calculable from the hybrid wavefunctions ( 110 ) and the 
values we have obtained are given in Table V. For the group-IV elements, the 
dielectric constnat is given directly in terms of / , V^ C and p / by Eqs . 

(51) and (55)- From Table II we calculate jb' = h 2 3 for £ ,= /, 3r for S; , 

- I,3i> for , and = /. 57 lor Sn. . For the binary compounds, the coupling 

matrix elements <£. Gf j i NO must be computed from Eq. (4-7). As remarked in 
Sec. II, the term ( (Gr( I* / ( E^“ E jE ) dominates the sum in Eq. (46) 

for the susceptibility % . We find 

|UV Eb )/ ' E T ) 1 1 ^61 I T > / <'$ I I ff > I * * a -»i (nit) 

in all cases. Finally, we scale our calculated values of %> for the group-IV 
elements in the manner of HC. Specifically, we replace in Eq. (51) by >v'f 
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and determine a for each row of the Periodic Table by fitting £ to experi- 

ment. We obtain Xu ^ 2,o2 for C » - - . » 2 for S; ,= 2. S i for ^ q , and \ h 
for S»u • For the binary compounds, we multiply our calculated by A a. A c , 
where A a. and- A £ are the appropriate values of A ^ for the anion and cation 
rows. Again as in the optical-absorption calculation and in contrast to HC, we 
do not use any direct experimental data for the binary compounds in determining 
the dielectric constant for these materials. In Fig. 2 we have plotted our- theo- 
retical dielectric constants versus the best available experimental values'*" 0 for 
the twenty materials under consideration here. The agreement with experiment 
is within 12$ for the binary compounds and, together with our results for £ , 

lends strong support to our method of determining parameters. 

We next consider the valence -electron density and the polarity of the bond. 
We have made a full evaluation of the single-bond electron density, as given by 
Eq. (60), in the cases of and Q-<xhs • These results are plotted in Figs, 

and i}, respectively. There is a broad peak in the electron density at the center 
of the bond in c.}e and a somewhat sharper peak near the Ac site in cm Ac . 

In both cases the bond is elongated perpendicular to the bond axis. The remainin 
peaks and valleys in the core regions result from the oscillations in the 
and jl atomic wavefunctions due to their orthogonality to the inner core states. 
We have also repeated these calculations in the Harrison limit using Eqs . ( 65 ) 
and (66) with o l c =; o( c and oi p = , the HC values of the covalency and 

the polarity. Interestingly, the results obtained are very close to those of 
the full calculations; in the bonding region the differences in both cases are 
on the order of 1$. This is consistent with the fact that we calculate o 4^ 
for <Ta, As as well as Q& (see Table V) and suggests that the electron- 
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density coefficients , U c and Ha.c are determined primarily "by the po- 
larity. As can he seen from Eqs . ( 65 } and (66), the latter is exactly true in 
the Harrison limit. 

To obtain the total electron density in the crystal, one must superimpose 
the individual single-bond densities. We have done this along a bond axis in 
both krtl and (quAs and the results are plotted in Figs. 5 and 6. The density 
near the center of the bond is increased only slightly by the overlap, but there 
is a significant increase in electron density in the back-bonding regions. For 

comparison we have also plotted in Figs. 5 and 6 the corresponding results of 

11 12 
Walter and Cohen obtained via the local-empirical-pseudopotential method. 

Their calculated densities are qualitatively similar to ours in the bonding and 

back-bonding regions, although clearly they find a higher electron density in 

the center of the bond than we do . 

An accurate experimental determination of the valence-electron density has 

been made in the case of Si , although not in either and A y . It has 
Ik 

been pointed out that the Walter and Cohen calculation for yields the 

magnitude of the central peak to within 7$ of experiment. In view of Fig. 5, 
it is of interest to consider what modifications in our Si 5 hybrids are required 
to in increase j’(F) in the center of the bond. From Eqs . (60), ( 65 ), (66) 
and (110), the electron density at the midpoint of the bond in a group-IV element 
is to a good approximation 

rui/z) = ~ [ M^/i) + 3 R r ( d /z)] 1 / 0 +s) . (115) 

Free-stom wave functions lead, to a value of 0.063 a.u. for f ( Iz ) in C ; , as 

13 

compared with the value of 0.102 a.u. found experimentally. However, by 
smoothly contracting the tails of R 5 ( r ) and beyond T =■ » one can 
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easily increase f ( 4(z ) to the required height 1 '*, as we have verified in a 
computer experiment. It seems clear, therefore, that valuable information 
about the shape of the hybrids can be extracted from an accurate knowledge 
of the electron density in the bonding region. This matter will be pursued 
elsewhere. 

Our defined polarity has been evaluated for each of the sixteen 

binary compounds considered here and these results, together with the carre- 

He 

sponding values of c/p , are given in Table V. For the non-nitrate hi , {fa. 
and I iv compounds, we find cAp approximately constant ( 0.51 to O.56), which 

He 

is roughly in accord with ol p ( 0.44 to 0 . 54 ). For the corresponding nitride 

,($• , HC 

compounds, however, we consistently calculate <Xp ^ -£• o(.p . In addition, 
^ „ c. r -on , ./ $ • >H< 


we 


find < °^p for S C and 5 W and c/p T > oLp for BP and BAs • 

These trends can be qualitatively understood in terms of our computed values 
of V 3 and . Generally speaking, cXp increases with increasing V i and 
decreases with increasing . We calculate V3 ^ and V3 ^ V4 > and 

consequently relatively small values of (Xp » for and all of the nitride 

compounds. For the remaining compounds, on the other hand, we obtain V3 > V, ^ d 

with V4 = Vj for BP and BAs an d V3 < V4 for the others. 

For comparison with oCp and o(p C , we have also listed in Table V our 
calculated values of 2?p > as defined by Eq. (69). Generally, there does not 

seem to be a simple correlation between 2 p and either (Xp or . We find 

(■p . ^ 

o(p for S; C and the nitrides, but < cXp for the rest. Interestingly, 


however, there is a qualitative correlation between (S c ~fa) and ol 


H £ 


In 


. , 1 He 

particular, the trend of decreasing cx.p with increasing mean atomic number in 
the B , Atf ■ > tfA. and I/n series (except for /U $i> ) are all reflected in 
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( £c.“ Scv'l . In this regard, note that ( measures the relative spatial 
extent of the cation and anion hybrids, while measures the relative weighting 

of these hybrids in the electron density. 

Let us now move on to consider the nuclear exchange and pseudodipolar 
coefficients. To evaluate Eqs . (85) and (86) for and Ppj , one needs values 
for the matrix element n and . Equation (80) for is easily 

evaluated in terms of the hybrids defined by Eq. (llO): 


)l W ■■ v " 


)in 


R Sti to) '[f? sn *(4)+ 3 (?pn"(:hj b h" 

[ R$h (d ) f ^ Rf>n (d ^ J 71 ~ n" A )i ' 

For free-atom wavefunctions >;> L +'3R^ n C't>j and consequently 

i." >>X. ■ Specifically, we find 

i /i,:: 1 < -• ' - 

U:' / i.;'„ 1 < °’ 01 

in all cases. The situation for is slightly more complicated. Only 
has a simple formula in terms of the hybrids: 


(116) 


(117) 


ft 


11 

U r' 


oo , A 

R* <r> <* 


(118) 


r 


To evaluate the overlap terms one must use the wavefunction-expansion techniques 

„ / ; 

discussed in the Appendix. Specific formulas for and are given there. 

Quantitatively, we find 



although the upper limits are only approached in the case of tJN . Typically 
the ratios are much smaller. For completeness we have included all four over- 
lap terms, and , in calculating P e and f' fJt . 

Values of P c and f 7 ^ have been determined both in the Harrison limit, 
using Y? and V 3 to evaluate Eqs . (93) and (9^)» and with the full theory, 
using the data listed in Tables III and IV to evaluate and E D . These 

results are presented in Table VI together with the best available experimental 
data. Note that the values of P e and in the Harrison limit are larger 
(except for I n /\j ) than the corresponding ones given in Table II of Paper I. 

This is due primarily to the appearance of the overlap matrix element $ in 
Eq. (93)> as was discussed in Sec. II. Also, as expected from Fig. 1, the effect 
of two-electron correlation is to enhance f e and in all cases. The magnitude 
of the enhancement, however, shows rather complicated trends depending both on 
bond length and polarity. In the group-IV elements, f -Ttf^for C , - 3, ? >o for 5/ 
and (it j and ~ 3.13 for Siv . 

We have not included any additional scaling factors in the theoretical numbers 
listed in Table VI . This does not appear to be important in the case of /’j,j , 
but clearly we overestimate the magnitudes of f' e , although the trends are correct 
We have made an approximate least-squares fit to the six non-zero experimental 
values of P e using a functional form 
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where A a and. A c are constants which depend only on row number in the 
Periodic Table and the fc are the full-theory values given in Table VI. We 
have thereby determined scaling factors A v , appropriate to each row: A,, - ' \ ■ / 
far the S; row, r t\ 1 7 for the (A row and r i\b c f for the 5» row. As with 
the scaling factor X K in the dielectric constant, A lv increases with increasing 
atomic number. Note, however, that An ^ ! while A „ I 

In Paper I we noted the following relationship between the dielectric 
constant and the nuclear exchange coefficient in the Harrison limit (without 
overlap ) : 

f - 1 = c r e Id , (121) 

where <2 is independent of and V 3 and depends only on intra-atomic para- 

meters. [See the discussion pertaining to Eq.. (7*0 of Paper I.] This motivated 
us to plot experimental values of - \ ) against the known experimental values 
of p , and we found a rather striking linear 

relationship in the series In P , J A/U and i„Sk . Since that time, experi- 
mental determinations of fg, and Pp.^ for A P > as well as new measurements on 

l6 

As', have been made by Cueman and Soest , allowing them to make a similar 
plot for the Qa. series. They have found a good linear relationship for that 
series too; both sets' of experimental data are shown in Pig. ?• We are now in 
a position to make a meaningful theoretical plot of (f •• ) ) vs. rv/Wi’ • 

Using our calculated values of £ from Fig. 3 and our full- theory values of f\ 
from Table VI, we have done this in Fig. 8 for all the III -V compounds . Ap- 
proximate linear trends can be seen in the J n andtfn. series, as well as the 
hi series, for the heaviest three compounds. To examine the origin of this 






behavior, we repeated the procedure in the Harrison limit. In this case, the 
linear behavior was improved slightly in the Ju series, worsened slightly in the 
fja series, and was destroyed in the Aj? series. Keeping otrly the most dominant 
terms in J" and , one has in li Hi ilarrison limit 


rYctc c - 



f mfi ( jvr (1-s‘f 

i (l-S 2 ) - ) 2 J 


( 122 ) 


/ r '^2 

where by construction V is the HC fitted value of fr* . Wow for 

the heavier <jo. and I>\. compounds , y 1 , £ , o(^ C and oP are rather 

constant. This implies that for a given cation series 




/v. C O » % t a „ £ 


(123) 


Evaluation of Eq. (123) gives 0.0155 forSk , 0.0163 for As , and 0.0311 for F . 

The higher value for P is consistent with Fig. 8, where the points for A? F > 

(sfn P and Jn P all lie above the straight line defined by the corresponding As and 
< l "pounds. 

A comparison of Table II of Paper I and Table VI shows that our full treat- 
ment of overlap and two-electron correlation has had only a modest impact on the 
important ratio Pp^ / P £ . Moreover, in no case has that ratio been increased 
significantly, which, as pointed out in Paper I, is necessary to simultaneously 
reconcile p ^ and | with experiment. It now seems clear that lj\j / f ' e is dominated 

by the ratio , which in turn is a direct property of the _S and £ 

wave functions which make up the hybrid. Interestingly, we have found in the case 
of S) that the latter ratio can be increased if the tails of the free-atom wave- 
functions are contracted in the manner necessary to account for the electron 
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density in the center of the bond. This matter bears further investigation in 
the case of binary compounds . 

As our final application, we turn to the cohesive energy. For the group 
IV elements, it is instructive to rewrite Eq. (104) as 

Ec^ ~ SO - 'u ^ Vi " * 4 - Edprr ~ ^ E'e* Epy~o , ( 124 ) 

where ~‘fE cc , t is the explicit contribution of two-electron correlation to the 
cohesive energy. Using Eq. (106) and comparing Eqs . (104) and (124), one has 
(with V L -0 ) 

E. X 2(h2S : )V/ | /(<-S J ) - Va H< - v 2 ) (125) 

- >. err 1 


As desired, Ecor r “ 0 for ~0 . Unfortunately, a complete evaluation 

of Eq. (124) requires the uncertain quantity A 4 a . Because an accurate 
evaluation of is quite difficult, we calculate only the remaining terms. 

Specifically, we replace in Eq. (124) with 



K + 


(126) 


and use our calculated values of V/' , V 2 , , $, and (\ to evaluate &ad 

audits components. The results are listed in Table VII together with the 
experimental values of the cohesive energy . Note that our calculated Eo/. 4 s 
less than experiment but positive in each case. This is consistent with small 

o 

negative values of n , as we expect theoretically. Also note that 
makes a large positive contribution and ~4 smaller, but non-negligible, 
negative contribution to Eat • We ma y contrast tnis with the calculation of E Co ^ 
made by HC. They, of course, neglected both E ^ and E Corf . Quantitatively, 
these omissions were partly compensated for by their use of smaller values of S 
(a constant value of 0.5) and a scaling factor of 0.8 multiplying Ep rp . Their 
prescription gives a rather fortuitous value of E a ,^- II.U-V for C , hut negative 
values of -2, 12 for E,‘ , “•VJ'f for Efe , and - o \ for . 
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We have determined empirical values of A £n for each group-IV element by 
fitting to the experimental cohesive energy: 


A <r 





) / ri . 


(m) 


As can be seen from Table VII, the magnitudes of thenf.,, are reasonable . The 
irregular variation from element to element is questionable and no doubt partly 
reflects the fact that our theoretical model is best for a large -band -gap 
material like C . In any case, we have used Eq. (lOj) to extend the calculation 
of to the binary compounds. We have assumed that a £ c + A £ ^ is a covalent 

quantity and have used the appropriate arithemetic mean of the group-IV values 
for each compound together with a direct calculation of the remaining terms in 
Eq. (100). The results are plotted against the known experimental cohesive 
energies in Fig. 9. The comparison is not as favorable as with E z and (- , 

but positive values are found in all cases except J>, ^ . 



IV. CONCLUSIONS 


We have presented here the full formal theory of the two-electron "bond- 
orbital model, extended its application, and extensively compared and contrasted 
it with Harrison's one-electron bond-orbital model. Both the analytic formulas 
given in Sec. II and the parameters listed in Tables I-V can be used to treat 
many additional physical properties, including those considered by Harrison 
and co-workers with the one-electron model. Our work here has shown that two- 
electron correlation effects are significant in the dielectric constant, the 
nuclear exchange and. pseudodipolar coefficients, and the cohesive energy, but 
perhaps not in the valence electron density, where the polarity and the shape 
of •$£' hybrid wavefunctions are the dominant quantities. In 6 - 1 , 

and , however, the effect of correlation is essentially multiplicative 

and is hidden to some extent if scaling factors, a • la HC, are used to achieve 
quantitative agreement with experiment. Nevertheless, it is encouraging that 
we have been successful here by actually fitting less parameters to experiment 
that did HC. An interesting and potentially important area for additional work 
is in the optimization of the s hybrid wavefunctions and in the impact this 

a rt _ 

has on calculated parameters such as V > V,, » a and P 

3 T H 
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APPENDIX 


We present here useful formulas for quantitatively evaluating the two- 
electron integrals U n and K and the other two-center integrals s 5 S c ,\ , 

and |V discussed in Secs. II and III. Central to the evaluation of 
1( ,, and K is the standard expansion 


t - r- 


= 'tu 




2JL* 


f ./ V- \y 

■r iw f) L(r) 


(Al) 


i>l 


where is the usual spherical harmonic and ( P?) is the lesser 

(greater) of p and p 7 . We choose the 2 axis to be in the bond direction 
such that the cation is located a distance 4 d 2 from the anion. Then only 
the 7>i - o component of Eq. (Al) (as well as any other similar expansion) 
will contribute to the integrals of interest, and one may work in terms of the 
Legendre polynomials Pj^ * 


Y» 0 ( r > = . 


V 4-tt 


(A2) 


In terms of the , the Sp 3 hybrid wave functions take the form 


t f ) “ “ L 


2 J*nr 


_ R sn (r) P P (^s9) •+ 3 0n JD P (Cose)] > 


(A3) 
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where p 0 { x 'i = j , P, (X ) = ^ » and Ot\ H and - - j , as above. 

Using Eqs . (Al) and (A3), the Coulomb potential arising from an j: - hybrid 


electron density can be written 

v„(?> = i t (r’> pr-vf h ( r ' > J r ' 

■ ! r - r I 

- 4 .[Zk(^ + 1 + 2©,2 s 'p(r^ p; M)4 z ff (r) (A4) 
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To evaluate K , one can expand V/ c <?Un spherical harmonics centered on the 


anion 


This type of expansion technique is discussed in detail by Lowdin 


17 


Again only the vn t= 0 component of the expansion is required and one can write 


tX> 


| ■ fe - o 

where j? is measured from the anion and 

k * 


(A9L) 


= fss, U. r'i + lfo. (4 , rl - if.,, . ^ r > ■+ 7 frr* 


4 s.'i w,r) 


(aio) 


with 


k 


pJ*'.k 


U, r) - 


2 k+ f 

r 




* i i 2 . 


f k ' l ~ 7 JF ) Ik ( rar (aii ) 


Id- r J 


Then one finds 


K 


- u 


f <|Vn V c ^(Pvi? 

1 ( j [Rs». ( r '] 2 + RW> Rf» <-rt (A12) 

+ [K f .(r)] 1 [lB=“^ rV+ 4 
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To compute the remaining integrals, one must use an expansion analagous to Eq. 
(A9) for the Sp 5 hybrid wave function: 


°n ^ 1 ^ n ^ ( t0 "' 


k *' 


£ j 


(AT3) 


where 


A» d, f) » ^ + fi (A rl 


(Al4) 


with 


o( 0 K (d,'r ) 


ZTt 


> / 2J+J_ 

4 r V 4-ir v q-rr 




id - r 1 




Formulas for $ , a and /p , in terms of the /\ ,, ( L •) follow immediately: 


S=:£\ r K s «. c-t Achlr) + Js Kf«Cr) Ac M.rljr'dr, (ai6) 


2 Sf j { .■f'.Rsix (r ^c^i r ^ + Kfo«r)l^(d,r) + ^AcC^nj (A17) 

o 


j rM 


and 


«*■ * 


m 



f- KsnOA A*. «.r> + S„ [ -f A,,. W,r) 


+ a,;. <-> 


I ‘V 


r 

<4 r 
r 
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i? *’ 

' Finally, a formula for 0. u may be derived, by making two expansions of the farm 


. x (whi ) 14 f A'- i-) | 

' 5 kX 

x | PJ*) p k ,-K) P k - •:* > ^ x , 


(A19) 


In this case an infinite, but rapidly convergent.; series is obtained. For each 
value of h,' there are one to three values of k which give finite contributions 
to the series. In our calculations we have accumulated terms up to = 3 and k -3. 

We omit results for the remaining Coulomb integrals f-j,, and J" , as these 
quantities subsequently drop out of analysis given in Sec. III. A formula for H h 
can be developed by using Eq. (Ak) and a single expansion of the form A13, In 
the case of ^ a double expansion as in Eq. (A19) is required and an infinite 
series is obtained. 
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TABLE I. Intra-atomic parameters for the group-III, -IV and -V elements 
obtained from the free-atom Hartree-Fock term values and wave- 
functions of Mann (Ref. 7)- Energies are given in eV, while 
^ n and are given in a.u. 








Element 

HF 

-£ 

n 


U 

n 


c 


B 

9.68 

1.26 

13.95 

17.7 

0.776 

Al 

6.96 

1.25 

9.27 

29.6 

1.09 

Ga 

7-14 

1.47 

9-37 

87.2 

2.89 

In 

6 . 56 

1.19 

8.43 

123 

4.46 

C 

13.14 

2.08 

17.62 

35.0 

1.66 

Si 

9.38 

1.80 

11.34 

48.1 

2.03 

Ge 

9.28 

1.96 

10.96 

120 

4.73 

Sn 

8.33 

1.57 

9.63 

160 

6.75 

N 

16.92 

3.09 

21.22 

60,5 

2.02 

P 

11.95 

2.42 

13.27 

7.13 

3.27 

As 

11.46 

2.48 

12.36 

157 

6.85 

Sb 

10.11 

1.97 

10.66 

200 

9.23 
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TABLE II. Intrabond parameters lor the group-IV elements and the group-IV 
and IIT-V binary compounds. The quantity Vp was obtained by 
fitting to experiment, while 8, K, and Vl ( were calcujuted 
using Hartree-Fock free-atom terra values and wave functions, us 
described in the text. Ail energies are given in eV; the bond 
length d is given in A. 


Material 

a 

S 

K 

v 2 

V 3 

v l, 

C 

1.54 

0.648 

14.70 

6.54 

0.0 

2.52 

Si 

2.35 

0.668 

9.83 

2.53 

0.0 

1.36 

Ge 

2.44 

0.659 

9.46 

2.42 

0.0 

1.33 

Sn 

2.80 

O.661 

8.37 

1.99 

0.0 

1.12 

SiC 

1.88 

0.627 

11.95 

4.07 

o.4o 

2.08 

BN 

1.57 

0.608 

l4.l8 

6.54 

2.27 

2.70 

BP 

1.97 

0.663 

11.65 

4.07 

1.74 

1.75 

BAs 

2.07 

0.656 

11.14 

3.98 

1.71 

1.77 

AIN 

1.89 

0.517 

11.27 

4.07 

2.33 

2.71 

A1P 

2.36 

0.633 

9.58 

2.53 

1.93 

1.4l 

AlAs 

2.43 

0.645 

9.34 

2.47 

1.93 

1.26 

AlSb 

2.66 

0.659 

8.68 

2.24 

1.63 

1.14 

GaN 

1.94 

0.513 

11.13 

3.98 

2.25 

2.83 

GaP 

2.36 

0.629 

9.58 

2.47 

1.84 

1.44 

GaAs 

2.45 

0.637 

9.28 

2.42 

1.83 

1.25 

GaSb 

2.65 

0.654 

8.70 

2.19 

1.54;' 

1.15 
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TABLE II 

(eont'd. ) 


InN 

2.15 

0.470 

InP 

2.5^ 

0 .602 

InAs 

£ .61 

0 .617 

InSb 

2.81 

0.644 


10.16 

3.61 

2.25 

2.99 

8.96 

2.24 

1.86 

i.)i8 

8.75 

2.19 

1.86 

1.33 

8.24 

1.99 

1.59 

1.12 
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Table III. Relative singlet and trip3,ct eigenvalues E® = K M~^ r 'a +r 'c~ K ^ 

for the twenty solids considered in the text, where = K^y and = 

0 

F*i Also listed are the theoretical predictions of the K? 

optical-absorption peak both from the present work, E^ = E^y - E^, und 
from the work of Harrison and Ciraci (Ref. 2), . All energies are 


in eV , 



Material 

F° 

e g 

E° 

e t 

E° 

E IV 

F° 

Theory 

P°_ F 0 p 0 HC 

E 2 E IY" E G E 2 

Experiment 8- 
E 2A E 2b 


C 

-2.32 

1*. 82 

9.87 

2h . 32 

12.2 b 

12.2 C 

12.2 

12.2 

Si 

-0.50 

1.19 

3.91 

9-98 

l*.l*l b 

1* . Uo c 

1* . 1*0 

1*. bo 

Ge 

-0.50 

1.15 

3.81 

9.51* 

l*.31 b 

l*.30 c • 

1*.3 

1* .3 

Sn 

-0.38 

0.89 

3.13 

7.88 

3.52 b 

3.52 



SiC 

-1.25 

2.33 

6.56 

15.63 

7.80 

7.98 

8.3 


BN 

-3.33 

1.78 

9.95 

2^.87 

13.3 

13.1* 



BP 

-1.81 

2.65 

6.0l* 

16.22 

7.85 

7.32 

6.9 


BAs 

-1.78 

2.55 

5.96 

15.89 

7.73 

7.2l* 



A£N 

-2.1*6 

2.23 

6.69 

17.08 

9.32 

9.07 



A£P 

-1.72 

1.32 

3.80 

11.28 

5.52 

1* -99 



AJIAs 

-1.80 

1. 39 

3.61* 

10.96 

5.M 

1*.85 



AJlSb 

-1.1*6 

1.20 

3.30 

9.83 

1* .76 

1* .68 

1*.25 

l* .6 

GoN 

-2.26 

2.06 

6.92 

16.90 

9.18 . 

9.26 



GaP 

-1.59 

1.22 

3.76 

11.02 

5.35 

5.11 

5.27 

5.74 

GoAs 

-1.71 

1.38 

... 3.60 

10.65 

5.30 

1*.93 

1**85 

5.33 

GaSb 

-1.36 

l.lfc 

3.26 

9.57 

l*.6l 

U.31 

l*.l 

1*.5 



6o 


Table III. (cont'd.) 


InN 

-2.13 

1.70 

6.6l 

15.96 

8.7!* 

8.5H 


inr 

-1.56 

1.02 

3. >18 

10. Ill 

5. Oh 

14 . 8)1 

1».8 5.1 

XuAs 

— 1. 6H 

1.07 

3.31 

10.12 

U .95 

I4.58 

1* .5 5.0 

InSb 

-1.39 

O. 9 O 

2.97 

9.01 

l».36 

l4.09 C 

I4.O8 


Reference 9 

b Fit to E^ C , as described in the text, 
°Fit to experiment 


vn 
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Table IV. Kxpansion coefficients entering Kq . (3l) for the singlet 
eigenstates M = VI (the ground state), M - IV and M = Y of the twenty 
Bemi conductors considered in the text. 

M=VI TV V 


Mate- 


rial 

a 4M 

a 5M 

a 6M 

n l»M 

ft 5M 

“6M 

n 4M 

a 5M 

a 6M 

C 

0.167 

0.167 

0.736 

0.928 

-0.920 . 

0 

-1.44 

-1.44 

2.33 

Si 

0.100 

0.100 

0.841 

0.950 

-0.950 

0 

-1.53 

-1.53 

2.47 

Ge 

0.103 

0.103 

0.838 

0.940 

-0.940 

0 

-1.49 

-1.49 

2.39 

Sn 

0.097 

0.097 

0.847 

0.942 

-0.942 

0 

-1.50 

-1.50 

2.41 

SiC 

0.198 

0.084 

0.781 

0.990 

-0.817 

-0.200 

-1.30 

-1.43 

2.15 

BN 

0.402 

0.020 

0.653 

1.11 

-0.574 

-0.642 

-1.0 6 

-1.48 

1.97 

BP 

0.465 

-O.O6I 

0.635 

1.23 

-0.512 

-0.865 

-1.21 

-1.71 

2.33 

BAs 

0.1*59 

-0.05& 

0.641 

1.22 

-0.508 

-0.855 

-1.17 

-1.68 

2.27 

AXN 

0.1* 1*8 

-0.015 

0.658 

1.10 

-0.357 

-0.814 

-0.680 

-1.32 

1.38 

AAP 

0.61*3 

-0.120 

0.494 

1.26 

-0.218 

-1.24 

-0.889 

-1.65 

1.92 

AZAs 

0.680 

-0.121 

0.451 

1.27 

-0.202 

-1.30 

-0.920 

-1.70 

2.00 

AJtSb 

0.659 

-0.131 

0.480 

1.30 

-0.238 

-1.28 

-1.00 

-1.75 

2.14 

GaN 

0.1*28 

-0.0 17 

0.678 

1.10 

-0.358 

-0.800 

-0.675 

-1.31 

1.36 

GaP 

0.619 

-0.121 

0.520 

1.25 

-0.231 

-1.21 

-0 . 886 

-1.63 

1.89 

GaAs 

0.655 

-0.113 

0.476 

1.26 

-0.224 

-1.25 

-0.904 

-1.66 

1.95 

GaSb 

0.632 

-0.128 

0.508 

1.29 

-0.259 

-1.24 

-0.998 

-1.72 

2.11 
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Table IV. (cont'd.) 


InN 

0.1*20 

-0.010 

0.695 

1.07 

InP 

0.626 

-0.118 

0.516 

1.22 

InAs' 

O.667 

-0.122 

0.1*71 

1.21* 

InSb 

0.671* 

-0.131* 

0 .1*66 

1.28 


-0.310 

-0.709 

-0.563 

-1.26 

1.16 

-0.179 

-1.22 

-0.761 

-1.55 

1.67 

-O.163 

-1.28 

-0.796 

-1.60 

1.76 

-o.ico 

-1.31 

-0.911 

-1.69 

1.98 
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Table V, Calculated parameters related to the dielectric constant and 


the polarization of 

the bond, 

as discussed in the 

text . 

The quantity 


Ct* 1 ^ is the 

P 

value of polarity 

assigned by Harrison 

and Ciraci (Re 

f. 2 ). 



Material 

6 

6. 

6 -6 

26 


G 

a 

HC 

a 

Y 


c 

a 

c a 

ca 

P 

P 

P 

1 

C 

0.279 

0.279 

0.0 

0.0 

0.0 

0.0 

0.0 

0.581 

Si 

0.265 

0.265 

0.0 

0.0 

0.0 

0.0 

0.0 

0.633 

Ge 

0.257 

0.257 

0.0 

0.0 

0.0 

0.0 

0.0 

0.646 

Sn 

0.259 

0.259 

0.0 

0.0 

0.0 

0.0 

0.0 

0.643 

SiC 

0.331 

0.229 

0.102 

0.098 

o.i4 

0.08 

0.39 

0.565 

BN 

0.334 

0.228 

0.116 

0.053 

0.24 

0.27 

0.41 

0.539 

BP 

0.274 

0.275 

-0 .001 

-0.040 

0.20 

0.35 

0.0 

0.602 

BAs 

0.261 

0.275 

-0.014 

-0.057 

0.19 

0.35 

0.0 

0.615 

AJIN 

0.390 

O.I89 

0.201 

0.110 

0.34 

0.34 

0.59 

0.492 

AJtP 

0.312 

0.229 

0.083 

0.054 

0.39 

0.53 

0.47 

0.593 

AMs 

0.303 

0.234 

0.069 

0.037 

0.40 

0.55 

0.44 

0.606 

A£Sb 

0.277 

0.251 

0.026 

o.oo 4 

0.36 

0.54 

0.54 

0.628 

GaN 

0.363 

0.184 

0.179 

0.114 

0.33 

0.32 

0.62 

0.528 

GaP 

0.299 

0.229 

0.070 

0.050 

0.38 

0.51 

0.52 

0.607 

GaAs 

0.288 

0.233 

0.055 

0.034 

0.38 

0.54 

0.50 

0.621 

GaSb 

0.266 

0.252 

o.oi4 

-0.001 

0.34 

0.52 

0.44 

0.638 

InN 

0.372 

0.166 

0.206 

0.156 

0. 36 

0.31 

0.64 

0 .523 

InP 

0.315 

0.213 

0.102 

0.096 

0.4l 

0.52 

0.58 

0.591 

InAs 

0.307 

0.218 

0.089 

0.080 

0.42 

0.56 

0.53 

0 . 6 o 4 

InSb 

0.285 

0.237 

0.048 

0.044 

0.39 

O.56 

0.51 

0 .625 



TABLE VI. 


% Vs 

. \ ' vv \ 


?-v 

^ -.-&* 

^ "y 4 
vr? V-' 




Nuclear exchange and pseudodipolar coefficients, T e and Fpa, for the twenty 
solids considered in the text in units of the direct dipole-dipole interaction 
coefficient F^. The theoretical results refer to our calculations done in both 
the Harrison limit and with the full two-electron formalism described in the text. 


i* 

t 



f. 

[•,: 

!:; 

9S: 

*■ 

r 


r e/ r dd r pd/ r dd r pd/ r e 

Theory Expt. Theory Expt. Theory Expt. 


Material 

Harrison 

Full 

Harrison 

Full 

Harrison 

Full 

C 

0 .0*4-3 

0.087 

0.0068 

0.0139 

0.160 

0.l60 

Si 

0.846 

2.79 

0.130 

0.427 

0.153 

0.153 

Gg 

6.00 

19.8 

0.860 

2.84 

0.143 

0.143 

Sn 

19. 7 

69.5 

3.04 

10.7 

0.154 

0.154 

SiC 

0.120 

0.435 

0.0187 

0.0681 

0.155 

0.157 

bn ■ , ; 

0.022 

0.055 

0.0034 

0.0076 

0.154 

0.137 

BP 

0.l6l 

0.203 

0.0254 

0.0386 

0.158 

0.190 

BAs 

0.404 

0.568 

0.0628 

0.103 

0.155 

0.180 

AIN 

0 .077 

0.273 

O.OO87 

0.0322 

0.113 

0.11c 

A1P 

0.412 

0.684 

O.O619 

0.102 

0.150 

0.150 

AlAs 

1.20 

1.5 1 * 

0.170 

0.214 

0.142 

0.139 

AlSb 

1.5b 

3.15 ~0.0 a 

0.233 

0.469 

0.149 

0.1-9 


o\ 

sr 


TABLE VI (cont'd. ) 


GaN 

0.185 

0.917 


0.0218 

0.107 


0.118 

0.117 


GaP 

0.913 

1.97 

o.i4 b 

o.i4o 

0.303 

0.20 b 

0.154 

0.154 

1.43 

GaAs 

2.69 

4.26 

0.73 b 

0.388 

0.623 

0.65 b 

0.144 

0.146 

0.89 

GaSb 

5 .07 

8.80 

1.89 a 

0.768 

1.36. 


0.151 

0.154 


InN 

0.277 

2.01 


0.0325 

0.246 


0.117 

0.122 


InP 

1.50 

4.00 

-0.55° 

0.228 

0.608 

~1.0 c 

0.152 

0.152 

-1.8 

InAs 

4.71 

8.72 

2 .06 a 

0.673 

1.24 


0 .143 

0.143 


InSb 

9.18 

15.6 

5.28 a 

1.38 

2.33 


O.15O 

0.149 


a R. K. 

Sundf ors , Phys . 

Rev. 185, 

458 (1969). 








Reference l6 

M. Engelsberg and R 1? Horberg, Phvc, Rev. B 5» 3395 (1972). 
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Tu'blo XII. Cohenive energy per atom pair and its components for the 

group - IV elements in eV. The quantities E° . and E . are our calcu- 

coll coll 

latcd values of the cohesive energy without and with the correlation 

energy E corr , respectively, as discussed in the text. The quantity 
fit/ 

Ac is the value of Ac necessary to fit E , to E CX ^'. 
n n coh coh 


Mate- 

rial 

8(l-s)v” C 

-4 E 

corr 

J ‘ E ol 

-E 

pro 

E° 

coh 

E T, 

coh 

£ expt 

coh 

a fit 

Ac 

n 

C 

17.18 

- 7.88 

21.44 

- 16 . 64 

22.0 

14.1 

14.7 ’ 

-0.075 

Si 

5 . 8 *i 

-3.84 

14.84 

-l4.4o 

6.28 

2.44 

9 . 28 ‘ 

-0.855 

Ge 

5.87 

-3.87 

14.24 

- 15.68 

4.43 

O .56 

7.74* 

- 0.898 

Sn 

4.77 

-3.23 

12.92 

- 12.56 

5.13 

1.90 

6.24 

-0.542 
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FIGURE CAPTIONS 


Fig. 1 


Fig. 2 


Fig. 3 


Fig. 4 


Fig. 5 


Correlation enhancement factor ; defined by Eq. (92) for 
the nuclear exchange and pseudodipolar coefficients of 
homopolar solids as a function of X ~ 2. V<^ /v? w ith \\ - n 
Two-electron bond -orbital -model values of the dielectric 
constant € .for sixteen binary compounds plotted against the 
best available experimental values (Ref. 10). No experimental 
data exists for BP » BAs or I^I\J and the estimates of Van 
Vechten (Ref. 10) have been used instead. 

Single-bond electron density for (qe as calculated from Eq. 
(60) using free-atom Hartree-Fock wave functions and the Cl ; 
given in Table IV. 

Single-bond electron density for iq-a. As as calculated from Eq. 
(60) using free-atom Hartree-Fock wave functions and the Cl ; ,■ 
given in Table IV. The As nucleus is on the right. 

Profile of the total valence electron density in along a 
bond axis, obtained by superimposing single-bond densities 
(Fig. 3) in the solid. The corresponding local-empirical- 
pseudopotential calculation of Walter and Cohen (Ref. 11) is 
shown for comparison (dashed line). The small arrows indicate 
the positions of the nuclei. 

Profile of the total valence electron density in (p a (\ $ along 
a bond axis, obtained by superimposing single-bond densities 
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Fig. 6 
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Fig. T 


Fig. 8 


Fig. 9 


(Fig. 4) in the solid. The corresponding local-empiric al- 
ps eudopotential calculation of Walter and Cohen (Bef. 11) is 
shown for comparison (dashed line). The small arrows indicate 
the positions of the nuclei with As on the right. 

Experimental values of 6 - 1 vs. experimental values of Te/I'y c\ ^ 
for the seven III-V compounds on which data is available. The 
linear trends in the In and 6r ft. series are to be compared 
against the theoretical predictions shown in Fig. 8. 

Theoretical values of € - I vs. theoretical values of /![<.) 4 4 
for fifteen III-V compounds, as calculated from the two-electron 
bond -orbital model. The corresponding experimental results, as 
known, are plotted in Fig. 7* 

Two-electron bond-orbital-model values of the cohesive energy 
Eco{ for sixteen binary compounds plotted against the known 
experimental values. References to the experimental data are 
given in the caption of Fig. 11 of Ref. 2. 
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